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A simple systematic rule, inspired by high-temperature series expansion (HTSE) results, is pro- 
posed for optimizing the expression for thermodynamic observables of ferromagnets exhibiting crit- 
ical behavior at T c . This "extended scaling" scheme leads to a protocol for the choice of scaling 
variables, r = (T - T c )/T or (T 2 - T 2 )/T 2 depending on the observable instead of (T - T c )/T c , and 
more importantly to temperature dependent non-critical prefactors for each observable. The rule 
corresponds to scaling of the leading of the reduced susceptibility above T c as Xc{T) ~ t~ j in agree- 
ment with standard practice with scaling variable r, and for the leading term of the second-moment 
correlation length as £*(T) ~ T~ 1//2 t~". For the specific heat in bipartite lattices the rule gives 
C* (T) ~ T~ 2 [(T 2 - T 2 )/T 2 ] ~ a . The latter two expressions are not standard. The scheme can allow 
for confluent and non-critical correction terms. A stringent test of the extended scaling is made 
through analyses of high precision numerical and HTSE data, or real data, on the three-dimensional 
canonical Ising, XY, and Heisenberg ferromagnets. For the susceptibility x(^) an d the correlation 
length £(T) of the three ferromagnets, their optimized expression, which consists of the leading term 
(respectively x*C^) an d £*(T)) and a quite limited number of confluent and non-critical correction 
terms, represents real data to surprisingly good approximations over the entire temperature range 
from T c to infinity. The temperature dependent prefactors introduced are of crucial importance 
not only in fixing the optimized expression at relatively high temperatures but also in determining 
appropriately the small amplitude correction terms. For the specific heat of the Ising ferromagnet, 
Cc (T) combined with two non-critical correction terms which are calculated with no free parameters 
once the correlation length critical parameters are known, reproduces real data nicely also over the 
whole temperature range. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

At a continuous transition, the expression F*(T) for 
the leading critical behavior of a thermodynamic observ- 
able F(T) has the well known form 



F C *(T) ~ (T-T c y 



PF 



(1) 



where T c and pp are the transition temperature and the 
critical exponent respectively. For the concrete analy- 
sis of numerical data, a normalization factor with non- 
critical behavior at T c must be introduced. The simplest 
and most traditional convention, which will be referred 
to below as T scaling, is to normalize each F*(T) by a 
temperature independent constant. For obvious reasons 
this constant is chosen to be Tg F for each observable; one 
then writes the normalized leading term as the familiar 
text-book expression: 



f:(t) =c f [{t-t c )/t c ]-p f =c F t-f p , 



(2) 



where t = (T — T c )/T c and Cf is the critical amplitude 
(see [l| for a detailed review) . An alternative and a priori 
equally valid choice is to write 



F*(/3) = C F [(p c -(3)/f3 c )-r 

PF 

C F 1 - & = C F r- pF 



C f [(T-T c )/T]-p f 



(3) 



where (3 is the inverse temperature 1/T and r = 1 — /3//3 c . 
Note that the temperature dependence of the normaliza- 
tion is now different for each observable. This "/3 scaling" 
form has become the standard normalization for theoret- 
ical work on the critical properties of ferromagnets and 
analogous systems, see for instance^" 3 -, although more 
complex normalizations have been used in special cases. 
At higher order, confluent and analytic correction terms 
(such as temperature independent constants) are intro- 
duced. Thus including the confluent correction terms, 
the critical behavior, F c ((3), is written in terms of the f3 
scaling as 



F c (0) = F*(/3)(l + aF r e - 
= C f t~ Pf (1 + a F T 6 



(4) 



where = voj with us being the [universal] confluent cor- 
rection exponent, and of is the confluent correction am- 
plitude. In the T scaling form, r in the above equa- 
tion is replaced by t. This critical scaling form is firmly 
established by field theory in the limit of temperatures 
very close to T^. Ratios of the of for different observ- 
ables are universal—. The exponent 9 is common in both 
scaling forms so long as < 1. However, no general 
argument seems to have been given which would show 
that either the T or the (3 scaling is optimal for all (or 
any) observables when a much wider temperature range 
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is considered. Recently we have proposed an extended 
scaling scheme for normalizing observables such that the 
leading critical expressions remain good approximations 
right up to the trivial fixed point at infinite temperature*. 
Our extended scaling scheme is based on a consideration 
of high-temperature series expansions (HTSE), and so is 
naturally formulated in terms of the (3 scaling. The most 
important ingredient of the scheme is the introduction of 
non-critical prefactors (3^ F in the normalizations, where 
each exponent <j>p is uniquely chosen such that the nor- 
malized F*(f3) tends to the correct asymptotic form in 
the limit T — ► oo. 

In the present work our aim is to further develop our 
extended scaling scheme to include explicitly the conflu- 
ent and analytical correction terms. We then validate 
our scheme by analyzing data for three canonical ferro- 
magnets: the 5=1/2 Ising, XY and Heisenberg mod- 
els on simple cubic lattices in dimension three. These 
models have been intensively studied over many years 
and their main critical parameters: T c , the critical expo- 
nents pf, 0, and certain critical amplitudes are known 
to high precision. Careful accounts of studies using dif- 
ferent complementary approaches are given for instance 
in Refs. Accurate simulation and HTSE results 

have been published in the form of tabulated data. The 
present analyses show that the appropriately normalized 
leading terms are good approximations over the entire 
temperature range, with small but identifiable correc- 
tions due to confluent and non-critical terms. We obtain 
estimates of non-universal critical parameters like crit- 
ical amplitudes Cp and confluent correction amplitudes 
a,F from the high precision numerical data. Our extended 
scaling analyses are in each case entirely consistent with 
field theoretical and HTSE estimates of the critical pa- 
rameters. 

An important result of the present analysis is to 
demonstrate that the prefactors (3^ F which have been 
introduced play a crucial role in extracting accurate val- 
ues of the critical exponents from simulation data even in 
a temperature range close to T c , such as t < 0.01. In the 
standard scalings without the prefactors the estimates of 
the leading critical term and of the confluent term from 
analyses of numerical data turn out to be modified to 
order ~ r (note t = r/(l — r)). 

The same approach based on the HTSE should be di- 
rectly applicable to a wide class of systems having the 
same intrinsic HTSE structure as the simple ferromag- 
nets. Extensions to more complex systems such as spin 
glasses are in principle straightforward^. 

The paper is organized as follows. In Scc.|TT]wc explain 
our extended scaling scheme for various thermodynamic 
observables, and discuss confluent corrections to scaling 
terms in our scheme. In Sec. IIIII we give methods of 
analysis for numerical data using our extended scaling 
scheme. We show how they work in practice for Ising, 
XY and Heisenberg fcrromagncts in Sec. IIVI IVl and IVI1 
respectively. In Sec. I VIII we make concluding remarks 
and discuss related problems. 



II. EXTENDED SCALING SCHEME 
A. Optimized expression for observables F{0) 
Let us suppose HTSE of an observable F(f3) is given 

by 

F(P) = a F ,o/3 0F (1 + a P ,i0 + a F , 2 2 + ■•■)■ (5) 

The most important ingredient of our extended scaling 
scheme is then to write F c (f3) as 

F c {(3) = R F {T)a Ffi (3^T-P^ (6) 

where 

R c f (t) =np(l + a F T e + ■■•), (7) 

with 1Z* F = Cp /(apfif3f F ). In particular, the leading con- 
tribution without the confluent correction is represented 
as 



F;{p)=TV F a Ffi P 



<t>F , 



(8) 



The idea here is to let F c (/3) not only represent the cor- 
rect power-law divergence t~ Pf with the critical ampli- 
tude Cp (and with certain confluent correction terms) 
at temperatures close to T c but also have an asymptotic 
form consistent with the HTSE in the high temperature 
limit. The observable F(/3) is then approximated as 

F{fl) ~ F&ifi) = F c {(3) + bp. /3^(l + b F . 1 (3+- ■ ■). (9) 

Here the second term represents the analytic (non- 
critical) correction term in the present scheme. Its co- 
efficients bp.o and bp,i are determined in such a way that 
Eq. ([9]), combined with Eqs. ([6]) and J7]), coincides with 
Eq. ([3]) termwise as a function of /?; for example, 

bp.o = a F , - p-^Cpil +a F + ■■■), (10) 



and a similar expression for bp^\. The above set of equa- 
tions with the minimum number of the confluent and 
analytic correction terms is an optimized expression we 
propose for the function F((3) which is analytic in the 
range < j3 < /3 C and is singular at /3 = C . An impor- 
tant quantity for analyzing our extended scaling scheme 
is Rf(t) defined by 



Rf(t) 



a F ,ot 



-PF 



(11) 



It is the ratio of the measured values of observable F(T) 
to its leading critical term including the (3^ F prefactor 
but without the critical amplitude Cf- Explicitly, in the 
vicinity of T c where F(f3) ~ F c ((3), it behaves as 

Rp{T)~R F {T)~np{l + a f T e + •••)• (12) 

The plot Rf(t) versus t 6 near r = thus becomes a 
straight line with intercept 1Z F and slope IZpaF, where 
the values of T c , pf and 9 are assumed to be known 
(cj)F and a,Ffi are given by HTSE analysis). In the limit 
j3 — > 0, on the other hand, it becomes Rf(t) = l + (op i — 
Pf/Pc)P + " " " ■ Between these limits the form of Rf(t) 
will depend on the entire collection of unspecified higher 
order corrections to scaling. 
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B. Susceptibility 

The "true" susceptibility, naturally defined through 
the magnetization response to an infinitesimal applied 
field, is given by the fluctuation-dissipation theorem as 

Xt{P) = Pjj'E(SiS j ) (13) 

ij 

The reduced susceptibility Xrod = Xt/P is (confusingly) 
almost always referred to in the literature as "the suscep- 
tibility" . For consistency we will follow this convention 
and write the reduced susceptibility as Xi but we will re- 
fer systematically in the text to "reduced susceptibility" . 

The HTSE for the reduced susceptibility x{P) m S = 
1/2 fcrromagnets is of the form with <f> x = and a x ,o = 1, 
or with abbreviation of a x .i = on, 

X(J3) = 1 + ai/3 + a 2 p 2 + a 3 p 3 + ■ ■ ■ . (14) 

Then the leading divergent expression, Eq. l[8]). is written 
as 

x*M = n x r-\ (15) 

with TZ* X = C x . The ratio R x (t) of Eq. (JTTJ) is reduced to 

R x (r)=x(P)/r-^ (=XM), (16) 

where X c(P) is Eq. © for x(T). Note that R x (0) = 
TL* X at T c , R x {t) = TZ* X (1 + a x r 6 + ■ ■ •) near T c , and 
R X (P) = 1 + (ai — 7//3 c )/5 + • • • near infinite temperature. 
If R X { T ) remains close to 1 over the whole temperature 
range (which is the case for the systems we consider as we 
will see below), the leading critical contribution without 
the correction terms, X* C (P) = 7?.*r _7 , is a good approxi- 
mation for the reduced susceptibility, X (P)- Furthermore, 
the small difference x(P) ~ xt (P) °f the Ising and XY fer- 
romagnets in the whole temperature range < p < (3 C 
turns out to be reproduced surprisingly well by our opti- 
mized expression, x° pt (P) °f Eq. ©, only with one con- 
fluent and two non-critical correction terms. 



C. Correlation length 

There are different alternative definitions for the cor- 
relation length, but any correlation length diverges at 
criticality as £(T) ~ (T — T c )~ v . The second moment 
correlation length £ sm is defined through the second mo- 
ment 

MP) = 5> 2 (SoS,) - 2d x (/3)£ sm (/3) 2 , (17) 

r 

with d the space dimensionality 3 -. From now on we will 
refer to £ sm (P) simply as £(/3). The HTSE results show 
that for N- vector S = 1/2 spins, the series for H2(f3) is of 
the form BiP+B2P 2 +B 3 p 3 + - ■ ■ and is well behaved with 



B\ = z/N, where z is the number of nearest neighbors. 
This yields fa = 1/2 and a €j0 = [z/2dN) 1 / 2 . We then 
reduce Eq. §8§ to 

^ = ^{mT T ~^ (18) 

where H\ = C ( /(zp c /2dN) 1 / 2 with C £ being the stan- 
dard critical amplitude in Eq. ([3|) for £. The non- 
standard normalization prefactor P 1 ' 2 for £*(/?) is our 
main result. The mean-field calculation^ of the corre- 
lation length through the fluctuation-dissipation theo- 
rem provides an example confirming the extended scal- 
ing form of Eq. (fl~8|) . See also the analysis of Fisher and 
BurfordiS, particularly their temperature dependent " ef- 
fective interaction range" parameter r±(T). 

The critically divergent part of £ (/3) with the confluent 
correction terms is represented by £ c (/3) and is written as 

UP) =ni(l + a,r e + ■■■) (j^j r~ v (19) 

The ratio R^if) becomes 

R*(r) = {m/{zP/2dN)^)/r- v . (20) 

Again, because of the confluent correction, it becomes 
R^t) ~ Pc /2 Kl(l + a^r 9 + ■■■) near T c and i? £ (/3) = 

1 + [{B 2 /2Bi) - (oi/2) - v/Pc]P H with a x being the 

coefficient in Eq. (|14[) near infinite temperature. 

D. Specific heat 

The usual analysis of the specific heat (defined as the 
derivative of the internal energy at fixed volume C (T) = 
dU(T)/dT\v) near criticality assumes the form 

C C (T) = C C [(T - T c )/T c ]- a + K (21) 

where a = 2 — vd and Cc is the critical amplitude of 
C(T), and it is standard practice to introduce a large 
non-critical (in fact temperature independent) contribu- 
tion K (sec e.g. Ref [U). 

While the series for the reduced susceptibility and the 
second moment /12 are polynomial functions of /3 with 
both odd and even terms, for bipartite (such as bec and 
simple cubic) lattices the HTSE expression for C(/3) con- 
sists of even powers of j3 onl y 12 ' 13 , and can be written 
as 

C(p) = c 2 p 2 + c 4 /3 4 + c 6 /? 6 + • • • . (22) 

One can carry through the same type of argument^ as 
in the case of fJ.2(T), except that as all the terms in the 
series are even in (/3//3 c ), the critical behavior must be 
re- written in terms of the scaling variable [1 — (P/P c ) 2 ] 
replacing [1 — (/3//? c )] in the equivalent expressions for 
the correlation length. Thus, with <pc = 2, one can write 
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the leading critical term, which corresponds to Eq. 
as 



c*(p) = p 2 n* c l- 




1 (T-T c 
2^2 



T 



(23) 

where H c = Cc2 a / ' P 2 - If the confluent correction terms 
are included, we obtain the expression corresponding to 
Eq. © as 



C c (/3) = C c *(/3) 



+ 2<>V P 2 



(24) 



where ac is the confluent correction amplitude. It is 
noted that, since the two critical amplitudes Cc and ac 
are introduced in the standard way (as represented by 
Eq. (HU) for Cc), the factors 2 a and l/2 e appear in the 
definition of 1Z C and in Eq. (|2"4")) . respectively. In fact 
there is a hyper-universal relationship linking this Cc to 



explicitly in the present paper, we note for reference that 
the extended scaling normalization modifies the FSS ex- 



pressions. The canonical FSS ansataii is 
F(T, L) ~ L pp / V F[L/£(T)}, 



(27) 



where F(x) is a universal scaling function. The fre- 
quently used FSS expression derived from Eq. I|27[). 



F(T, L) ~ L PF ' V {T - Tc)\ 



(28) 



contains the implicit assumption of T scaling for the cor- 
relation length. It is thus only appropriate if restricted 
to a very narrow range of temperature around T c . With 
the extended scaling and the finite size correlation length 
£(L,(3), the FSS ansatz can be rewritten^ 



F(L,0)~ p** 



L 

W 2 



Pf/v 



T 



P 1 / 2 



1 - 4r 



(29) 



-'hyper 



(aC c ) 1/d C ( = 



faK* c 



V 2 a 



1/d 



(2d) 



1/2 



(25) 

where Chyper is a constant whose value is known rather 
accurately-^. Equation (f2"3"|) is not standard, but it can 
be seen to tend to the appropriate limit, C*(T) ~ (T — 
T C )~ Q , as T approaches T c . 

In practice H* c is much larger than unity (as will be 
seen later in the case of the 3d Ising model) which is the 
reason for the large non-critical contribution to C(/3). 
The non-critical contribution is in fact not a parameter 
to be adjusted freely, but it has to be determined through 
the high temperature limit of an equation which corre- 
sponds to Eq. ([9]). Ignoring the confluent correction so 
as to clarify the discussion, we know the exact high tem- 
perature limits for C c (/3)(~ C**(/?)) from Eq. and for 
C(J3) by Eq. fl52J). Then C(0) truncated to two leading 
non-critical correction terms is explicitly written as 



COS) -C c (/3)+/C 2/ 3 2 +/C 4 /3 4 , 



(26) 



where the non-critical parameters K-i are given by /C2 = 
c 2 — Tic an( l ^4 = C4 — all c I (3 2 . The coefficients C2 and 
C4 are known from HTSE. So if v, (3 C and i?| have been 
measured independently, we can evaluate all the param- 
eters which one needs to fix the functional form of C(/3), 
such as a = 2 — vd and H* c determined through Eq. (f25|) . 
This we assume to be C opt (/3), an optimized expression 
for C(/3), in the whole (3 range < j3 < f3 c . The thus cal- 
culated curve C opt (/3) can be tested by comparing with 
simulation and HTSE data. 



E. Finite size scaling 

Though we will discuss thermodynamic limit behavior 
only and will not analyze finite-size-scaling (FSS) data 



F{L,0) 



L 



T 



L 



(30) 



where the scaling functions behave as !F(x) ~ x PF and 
F(x) ~ x~ pp l v at x 3> 1. For the reduced susceptibility 



with 



0, the FSS form is written as 



L 



T^F-y 



F, 



[3 1 / 2 



l/u 



1/u 



Pc 



(31) 



where F x (x) ~ x 1 and F x (x) ~ constant, at x 3> 1. In 
a similar manner, the FSS form for the correlation length 
£(L,P), for which p^ = v and <j>^ = 1/2, is written as 



£(L,P)~LF £ 



L 

pv 2 



l/u 



Pc 



(32) 



where F^(x) ~ x~ v at 1 > 1. While the extended 
FSS scheme for the susceptibility is modified from the 
standard one only by the /3-prefactor in the argument of 
F x (x), the scaling plot is significantly improved for 2d 
Ising ferromagnetic and 3d Ising spin glass models^. 



III. ANALYSES USING EXTENDED SCALING 

In order to make a stringent test of the extended scal- 
ing scheme, we study the three canonical ferromagnets: 
Ising, XY and Heisenberg, on three dimensional simple 
cubic lattices. High precision numerical data have been 
obtained for each of these systems for the temperature 
domain ranging from close to T c to about 1.1T C and the 
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authors have generously published their data in tabulated 
for m 16 i 17 ' 18 . The data have been taken on systems large 
enough for the data points to be representative of the 
thermodynamic limit. Long HTSEs have also been pub- 
lished for x and A*2 and for C for all three system o 12 ' 19 , 
and relatively longer series for the free-energy and the 
specific heat have been calculated for the Ising modeli^; 
these series can be used to calculate x(T), £(T) and C(T) 
explicitly for the region T well above T c . Below we call 
these HTSE and MC data as the real data. Thanks to a 
combination of results from field theory and HTSE the 
values of the critical temperatures, the critical exponents 
and the critical amplitudes are known to a high degree 
of accuracy, and the confluent correction exponents are 
also well known. The [non-universal] confluent correction 
amplitudes are small for these three systems and the es- 
timates are much less accurate (see Butera and Comii 9 - 
for a detailed account). 

In each case we will plot the ratios R x (t) and R^{t) 
respectively defined by Eqs. fTB")) and ([20)1 with respect to 
t s , where we have used z = 2d for simple cubic lattices. 
The plots near r = give us the critical and the confluent 
amplitudes as explained at the end of Sec. Ill Al The 
ratios arc defined in a way that they approach unity at 
infinite temperature. We will see that the two ratios 
defined in the extended scaling scheme are in fact close 
to unity (within several hundredth deviation from unity 
at most) in the whole range of r. 

In addition, a simple scaling relation links the observ- 
ablcs xi§) and £(/?) through xW ~ £ 2 ~ r '(/3) to leading 
order. This equation has the advantage that it can in 
principle be used to determine the exponent r\ directly 
from a log-log plot of x{P) against £(/3) near f3 c without 
any explicit knowledge of (3 C . For the extended scaling 
scheme {(3 scaling with the (3^ F factors) , the relation can 
be rewritten to leading order 



X(P) 



2-7? 



We will analyze the ratio defined by 

x(/3) 



R(r) = 



2-/3 ' 



(33) 



(34) 



Including the leading confluent correction factors, it be- 
haves near T r as 



R(t)~K*(1+Bt + •••), 



(35) 



where K* = ^*/(^|) 2 -" = C X /(Q//3 C 1/2 ) 2 -" and B = 
a x — (2 — ri)a^. We also note that the ratios a^/a x are 
universal and are known to be about 0.7^2. This means 
that B 0.40a Y . 



IV. 3D SIMPLE CUBIC ISING FERROMAGNET 

For the 3d simple cubic Ising case N = 1, together 
with the high precision MC data at temperatures close 



n: io 




FIG. 1: An extended scaling plot of \ against in the 3d 

Ising ferromagnet. The filled squares represent the high pre- 
cision MC data by Kim et a^ and the filled circles numerical 
estimates from the HTSE of Butera and Comii^. Monte Carlo 
data with L = 16, 32 and 48 by ourselves are also shown. The 
straight line has a slope of 2 — r\ with r\ — 0.037(1). In this 
and the following figures our MC data are finite size limited 
for T close to T c , particularly in the case of L = 16. 
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FIG. 2: A conventional scaling plot of \ against £ in the 3d 
Ising ferromagnet. The data are the same as in Fig. [TJ 



to T c by Kim et aU^ and the HTSE estimates at relatively 
high temperatures by Butera and Comii^, our own MC 
data are also used in order to interpolate them and to see 
overall temperature dependences of x>£ and C. In our 
simulation we used the exchange MC method in combina- 
tion with 64 bit multi-spin coding technique for making 
equilibration fast. The 64 different temperatures simu- 
lated are distributed in the range of 4.0 < T/J < 15.0. 
The amount of total MC steps for L = 48 is 2.4 x 10 5 
and the last 8 x 10 4 MC steps are used for taking thermal 
averages. 

Figure [1] shows the parameter free log-log plot in the 
extended scaling form of the reduced susceptibility x 
against £/y/]3 data. Without allowing for corrections, 
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the slope of the line fitted to the data points (ignoring 
our MC data when they are polluted by finite-size ef- 
fects) gives a first estimate rj ~ 0.037. Figure [2] is the 
equivalent standard (T or (3 scaling) log-log plot of \ 
against £ with the slope fixed to the one obtained from 
Fig. [1] It can be seen that in the standard scaling form 
the linear relationship breaks down rather quickly while 
in the extended scaling form with the same input data, 
the linearity persists to a good approximation up to an 
infinite temperature and down to temperatures near T c 
until limited by finite-size effects. 

We examine the leading correction of the extended 
scaling formula given by Eq. (|35p . To higher pre- 
cision, Fig. [3] shows a plot of x/(S,/VP) 2 ~ v against 
t\ assuming f3 c = 0.2216544, rj = 0.0368 and 9 = 
0.50 4 7 ' 20 . The line is obtained by fitting the data points 
at t 6 < 0.4 to Eq. i[3"5"]). The intercept at r = 0, 
1Z* = 0.971(4), is in good agreement with the value 
C x /(Cf/(/3 c ) 1 / 2 ) 2 -'' = 0.9767(20) assuming the critical 
amplitudes from HTSE 19 . From the initial slope, B = 
a x — (2 — T])a^ = 0.086(11), which we will comment on 
below. 



at small r, we obtain TZ* X = 1.132(6), K\ = 1.074(3), 
a x = -0.138(23) and a £ = -0.109(20). The TZ* F values 
are in excellent agreement with the HTSE estimates^, 
K* x = C x = 1.111(1) and n\ = C ( /f3c /2 = 1.0677(7). The 
a f values are in qualitative agreement with the HTSE es- 
timates a x = -0.10(3) and a fi = -0.12(3)12. 
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FIG. 3: A plot of x/(UVP) 2 ~ v against t° in the 3d Ising fer- 
romagnet. The values of the critical parameters are assumed 
as described in the text. The straight line is the fit to Eq. (|35[) 
with TZ* = 0.971 and B = 0.086. 

Figures IH and [S] show the ratios R x (t) and R^(t) of 
Eqs. lfiTj|) and (f2TJ)) . respectively. The numerical data are 
taken from Kim et ali£, and the higher temperature val- 
ues are calculated using the tabulated series of Butera 
and Comii 2 -. The HTSE terms were simply summed, 
and the points quoted correspond to the temperature 
range where the contributions from further terms can 
be considered negligible on the scale of the plots. By 
using appropriate extrapolation techniques, like differen- 
tial approximations, the range over which the published 
HTSE datapi 2 - could be used to evaluate the temperature 
dependence of the observables to high precision could 
be considerably extended. The assumed critical param- 
eters are 7 = 1.2372, v = 0.6302 and 6 = 0.504 ^. 
From the initial intercepts and slopes of the fitted line 



FIG. 4: A plot of the ratio R x ( r ) against r in the 3d Ising 
ferromagnet. The straight line represents a fitting to R X ( T ) = 
TZ* X (l + a x T e ) with TZ X = 1.132(6) and a x = -0.138(23), while 
the curve does R x (r) calculated from x° pt (P) of Eq. © . In 
the inset, the T scaling ratio R x {t) against t e is shown. 

An overall conclusion on the extended scaling analy- 
sis of the 3c? simple cubic Ising data, which will be con- 
firmed by the analyses of the two other systems as well, 
is that this form of scaling is entirely consistent with the 
high precision values of critical parameters from exten- 
sive HTSE and field theoretical (FT) work. It is remark- 
able that over the entire temperature range from T c to 
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FIG. 5: A plot of the ratio R^(r) against r in the 3d Ising 
ferromagnet. The broken straight line represents a fitting the 
data by Kim et al. to R((t) = TZ^(l + a^r e ) with 7?.^ = 
1.074(3) and a ( = -0.109(20), while the curve does R^r) 
calculated from £° pt (/3) of Eq. ©. 
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infinity, the maximum deviations from the leading criti- 
cal expressions of Eqs. (fT5j) and (|18[) are of the order of 
a few percent. Let us go into further discussions about 
the R x (t) behavior. In the inset of Fig. 2J we show the 
corresponding T scaling ratio R x {t)=\{T) / 1 _1 plotted 
against t e . The latter is calculated using the same val- 
ues of the critical parameters T c , 6, 7 and C x as those for 
R x (t), and so by construction in the low i, r limit the 
intercepts and slopes of both ratios must coincide. It can 
be seen that in fact the T scaling curve only approaches 
the (3 scaling curve closely in the range of t, t extremely 
close to zero. This result for x with <f> = strongly 
suggests the superiority of the (3 scaling, and hence our 
extended scaling, over the T scaling. 

The full curve in the main frame of Fig. 0] is the opti- 
mized expression, i?° pt (r), which is evaluated through 
X° pt (r), with one confluent correction term discussed 
above and the two non-critical terms. The first term of 
the latter is a constant, 6 Xi o = 1 — C x , which yields sim- 
ply i?° pt = 1 at r = 1, or at an infinite temperature. Its 
second term which is also calculated via the parame- 
ters already fixed, specifies the slope of i?° pt (r) at r = 1 . 
By taking into account only these three correction terms 
to the leading critical term, \t (P) of Eq. (fT5"|) , we obtain 
i?° pt (t) which reproduces surprisingly well the real data 
in the whole temperature range < r < 1. Notice that 
t 6 = 0.6 corresponds to T ~ 1.57T C . This result also in- 
dicates the superiority of our extended scaling with the j3 
scaling: xt (/?) not only represents the critical behavior of 
x(/3) close to T c but also x(P) m the whole temperature 
range up to infinity. In this context we note again that 
= 1 for the "true" susceptibility and that the reduced 
susceptibility %(/?) is derived through our extended scal- 
ing scheme. We also note that the similarity between the 
R x (t) plot in Fig. 3] and the R^(t) plot in Fig. [5] over the 
entire range of temperature is striking. 

Lastly, Fig.[Ushows C{f3) / (3 2 as a function of l-/3 2 //3 c 2 ■ 
The data points are calculated from the HTSE of Arisue 
and Fujiwara which extends to powers up to 2?i = 46^, 
MC energy data at L = 128 and 96ii, and our numerical 
simulations for different sizes up to L = 48. We examine 
the extended scaling with non-critical contributions to 
C(P) given by Eq. ((26)) . By using the hyper-universal re- 
lation with the value of Ch ype r equal to 0.2664(1)2 for the 
3d Ising model and our £ analysis, we obtain TV C ~ 29.4. 
Then the non-critical parameters Yii and /C4 are deter- 
mined by C2 and C4 of HTSE and with putting ac — 0. 
The solid curve represents the no-free parameter plot 
of Eq. (|2"6"|) with the a, 1Z C , C2, C4 values cited or esti- 
mated above. The agreement over the whole tempera- 
ture range is very satisfactory; the non-critical correction 
is so strong that the bare leading power law is a poor ap- 
proximation until very much closer to T c than the range 
covered by the figure. We consider this result as an indi- 
cation that the extended scaling scheme combined with 
the optimized introduction of correction terms is an effec- 
tive method for analyzing critically-divergent quantities 
in general. 




0.01 0.1 1 

i-(/m) 2 

FIG. 6: A plot of C(f3)/P 2 against 1 - P 2 /Pt- The filled tri- 
angles represent the numerical estimates by HTSE of Ref. [3| 
and the filled squares the MC data of Ref. [lU • The solid line 
represents the expression of Eq. (|26[) with TZq = 29.4 and 
ac = 0.1. The straight broken line is the bare leading critical 
power law as (1 - (/3//3 c ) 2 )~ a . 




FIG. 7: An extended scaling plot of xW) against £/-v/ ' P/2 
in the 3d XY ferromagnet. The squares represent the 
high precision MC data by Gottlob and Hasenbuschi£ and 
Hasenbuscb^i, and the filled circles the numerical estimates 
by HTSE of Butera and Comii^. The fitted straight line has 
a slope of 2 — r\ with r\ — 0.036. 



V. 3d XY SIMPLE CUBIC FERROMAGNET 

The same analysis has been carried out for the 3d XY 
model (N — 2). High precision numerical data were 
published by Gottlob and Hascnbuschii, and are sup- 
plemented here by unpublished data kindly provided by 
M. Hasebusch^i. The higher temperature data are cal- 
culated using the tabulated series of Butera and Comii^. 
The critical point is (3 C = 0.4541652(5) and the exponents 
77, 6, 7 and v arc close to 0.0381, 0.53, 1.3178 and 0.6717, 
respectivel y 7 ' 8 ' 19 . Figure [7] shows the x ~ / (/3/2) 1 / 2 log- 
log plot. The leading scaling scheme works well up to 
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FIG. 8: A plot of x/H/xfW/^f^) against t° in the 3d XY 
ferromagnet. The critical parameters are assumed as f3 c — 
0.4541652, T) = 0.0381 and 6 = 0.53. The solid line shows a 
fitting line to Eq. d35j with TV = 0.9995(2) and B = 0.047(1). 



very high temperatures, as in the Ising case. The slope 
in Fig. [7] gives us the value of 77 which is in agreement 
with the previously rep orted values^. Figure [5] shows 
the plot of x{ T )/[ii. T )l \fWFA 2 ^ 71 against t 6 assuming the 
central values for the exponents 77 and 9 as mentioned 
above. Figures 151 and [TTjl show R x {t) and R^(t) respec- 
tively against r e . From the r = intercept and the initial 
slope one can estimate K* x = 1.0471(4), K\ = 1.0238(3), 
a x = -0.093(3) and a 4 = -0.073(2). These are all rea- 
sonably close to the quite independent HTSE values^ 
K* x = 1.014(1), TZ\ = 1.0102(6), a x = -0.04(2) and 
ci(. = —0.07(3), but are probably more reliable as they 
are consistent with the independent FT estimate of the 
universal ratio a^/a x ~ 0.65, see comments in Ref. fl9| . 
Also, the values of R* and B in Eq. (|3"5"|) calculated from 
thus obtained set of the parameters reproduce well the 
data as shown in Fig. [5] This agreement again validates 
the extended scaling protocol and demonstrates that a 
combination of information from FT, HTSE, and simu- 
lations analyzed using this protocol can lead to consistent 
high precision critical parameter measurements. 

For comparison, we plot the standard T scaling ratio 
R x (t) introduced in Sec. IIVI also in Fig.[5J Its coincidence 
with R x (t) will only hold for t <C 0.01. As is the case for 
the Ising system, the slope of R x (t) is opposite to that 
of R x (t) and the magnitude of R x {t) — R x {0) is much 
larger than the corresponding magnitude of the extended 
ratio already at t s ,T e ~ 0.2, or i,r ~ 0.04. In Fig . [TO] 
we also show the T scaling R^t) = (£(T)/ yJ(3 c /2)/t~ u 
and the i? c (r) = (£ (T) / i/ ( 9 c /2)/r _I/ by /^-scaling. The 
true leading term plus confluent correction holds with 
the extended scaling form, R^(t) of Eq. (f20j) with N = 2 
up to t ~ 0.1 while with the other forms of scaling the 
correct limit will hold only for t -C 0.01. In particular, 
the comparison of (3 scaling R^ (r) with extended scaling 
R$(t) demonstrates the importance of the /3 1 / 2 prcfac- 



FIG. 9: A plot of the ratio R x (t) against t in the 3d XY 
ferromagnet. The value of 7 is assumed to be 7 = 1.3178^. 
The filled marks are obtained by the extended scaling of 
HTSE i2 and the MC data^i. The broken line represents a 
fitting to R x {t) = 7l* x {l + a x T e ) with Tl* x = 1.0471(3) and 
a x — —0.093(3), while the curve does R X ( T ) calculated from 
X opt (/3) of Eq. ©. The standard T-scaling ratio R x (t) as a 
function of t e is also plotted by the open circles. 

tor in Eq. (|2T)|) of the extended scaling scheme. These 
results imply that even close to T c the extended scaling 
is a considerable improvement over the standard scaling 
analysis for estimating critical parameters including the 
correction terms. 

The curve in Fig. [^represents our optimized estimates 
i?° pt (r) up to the second order of non-critical corrections. 
It reproduces about 5 percents change in i? x (r), from 
about 1.05 at r = to 1 at r = 1, to a very good approx- 
imation. The corresponding relative change in R(-(t) is 
only less than 2 percents as seen in Fig. [IJJ To reproduce 
this change by i?^ pt (r) to an approximation as good as 
i?° pt (r) in Fig. [9l however, more than third order non- 
critical correction terms are required. 



VI. 3d HEISENBERG SIMPLE CUBIC 
FERROMAGNET 

The same analysis has been carried out for the 3d 
Heisenberg model (N = 3). High precision numerical 
data were published by Holm and Jankei^, and are sup- 
plemented here by higher temperature data calculated 
using the tabulated series of Butera and Comii^.. The 
critical point is f3 c = 0.69305(4) and the exponents 77 and 
9 are close to 0.036 and 0.55^^. A recent exponent set^ 
gives 7 = 1.3960(9), v = 0.7112(5) and 77 = 0.0375(5). 

Figure HD shows the X {T) - (,{T)/^/Jj3 log-log 
plot, which gives an estimate of 77 consistent with 
that of Ref. [22j. Figure [12] shows the plot of 
x(T)/[£(T) / yJ '(3/3] 2 ~ v against t s assuming the exponent 
values as 77 = 0.0375 and 9 = 0.55. From this plot it 
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FIG. 10: A plot of the ratio R^(t) against r" in the 3d XY 
ferromagnet. The value of v is assumed to be v — 0.6717. 
The line represents a fitting to R^(r) = 7?.|(1 + a^r 6 ) with 
Til = 1.0238(3) and a 5 = -0.073(2). The ratios from T- 
scaling R^(t) and /3-scaling R((t) are shown by open circles 
and open triangles, respectively. 
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FIG. 11: An extended scaling plot of x{P) against £/ \/ /3/3 in 
the 3d Heisenberg ferromagnet. The filled squares represent 
the high precision MC data by Holm and Janke^ and the 
filled circles the numerical estimates from the HTSE of Butera 
and Comii^. The straight line has a slope of 2 — r\ and the 
best fit gives T] = 0.0379(4). 



appears that the initial slope is very small, correspond- 
ing to almost zero values for a x and ag. Figures [T3l 
and [T4l show respectively R x (t) and Rf (r) against t 6 , 
assuming the values of 7 and v in Ref. [22|]. The MC 
and HTSE points may not appear to connect smoothly 
in these figures, because the manner in which the plots 
are presented enhances small deviations from the lead- 
ing term form. However, the change in the values of 
both R x (t) in Fig. QJ and R^(t) in Fig. M are lim- 
ited to within a few percent of their absolute magnitude 
in a whole range of r as is the case for the other two 
ferromagnets studied. From the straight line fit of the 
MC data at small t 9 , one can estimate 1Z* = 0.952(2), 
Kl = 0.967(2), a x = -0.04(1) and = -0.03(2). In 
this case the parameters are slightly less consistent with 
the HTSE estimates^, K* x = 0.9030(8), K\ = 0.9447(5), 
a x = 0.06(3) and = 0.003(6), but it should be noted 
that the estimates for these [non-universal] parameters 
depend very sensitively on the precise values taken for 
the critical exponents. We certainly need more precise 
data near T c to fix the values of these critical parameters 
for the Heisenberg ferromagnet. 



VII. CONCLUSION 

We have outlined a systematic rule for the scaling 
and normalization of thermodynamic observables having 
critical behavior at continuous phase transitions. This 
"extended scaling" rule corresponds for ferromagnets to 
scaling of the leading term of the reduced susceptibility 
above T c as Xc{T) = 7?.*t~ 7 in agreement with stan- 
dard practice, for the leading term of the second moment 



correlation length as f c (T) = TZW 1/2 t- 



with Tl\ = 
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FIG. 12: A plot of xAC/V/VS) 2 "") against t 6 in the 3d 
Heisenberg ferromagnet. The critical parameters are assumed 
as f3 c = 0.69305, 77 = 0.0375 and 9 = 0.55. 



Q/ ' {z(3 c /2dN) 1 / 2 and for the leading term of the specific 
heat in bipartite lattices C(T) = {3 2 1Zl 



1-Nr- 



A 

with n* c = C c 2 a //3 2 plus strong non -critical correc- 
tion terms which we explicitly evaluate by linking to the 
HTSE. 

Analyses are made of high precision numerical data 
on three canonical ferromagnets using these expressions 
allowing for confluent scaling correction terms, plus non- 
critical corrections for the specific heat. Near T c the re- 
sults are entirely consistent with the critical parameter 
sets (including the confluent corrections) which have been 
obtained independently using sophisticated FT, HTSE 
and simulation techniques^ i 19 ' 20 i 22 . 
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FIG. 13: A plot of the ratio R X ( T ) against r in the 3d 
Heisenberg ferromagnet. The straight line represents a fitting 
to R x (t) = ft* (1 + a x r e ) with ft* = 0.952 and a x = -0.04, 
while the curve does the optimized form using up to the sec- 
ond non-critical correction terms. 

1 , , , , , ^ 

0.99 - • 
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FIG. 14: A plot of the ratio R^(r) against t° in the 3d 
Heisenberg ferromagnet. The line represents a fitting to 
R s (t) = ft|(l + a s r") with ft* = 0.967 and a $ = -0.03. 



The most important result found in the present work 
is that, for x(T) and £(T) the leading critical expres- 
sions with the extended scaling normalizations F*(j3) of 
Eq. ([5]) agree to a very good approximation with the true 
F{(3) up to infinite temperature. To demonstrate this 
fact more in details we have introduced the ratio Rf(t) 
defined by Eq. (jlip . For \ of the Ising ferromagnet, for 
example, it is equal to the critical amplitude C x at T c 
(r = 0) and to unity at infinite temperature (r — 1) by 
definition. R x (r) evaluated from the true data are repre- 
sented by the data points in Fig. 31 while R x (t) evaluated 
through the leading expression x*((3) is independent of r 
and equal to C x . The difference between the two is, how- 
ever, at most 13 percent in this case. The corresponding 
differences for R x (t)'s of the two other ferromagnets as 
well as for R^ (r) 's of the all three ferromagnets are less 



than several percent. This is our first result mentioned 
just above. 

We have next demonstrated that our extended scaling 
scheme, in terms of the (3 scaling and with the tempera- 
ture dependent prefactor j3^ F , is of crucial importance in 
precisely extracting the small amplitude ap of the lead- 
ing confluent correction term. The result is represented 
by the solid line in Fig.[4]as well as those in Figs. El GO and 
[TU1 In addition, we have also checked that the optimized 
expression F opt (/3) of Eq. © , consisting of F c * (/3) and one 
confluent and two non-critical correction, yield i?^ pt (r) 
which reproduces the true Rf(t) surprisingly well as 
shown the curves in Figs. Q] and [51 though more than 
third non-critical correction terms would be required for 
equally good agreement in other observablcs. 

The large non-critical terms in the specific heat C(T) 
are also incorporated explicitly within our extended scal- 
ing scheme with no further adjustable input parameters. 
For the Ising ferromagnet on the simple cubic lattice 
C(T) is calculated to a good approximation over the en- 
tire temperature range (sec Eq. (|26|1 ). Although the non- 
critical correction terms are large for C (T) , the principle 
of the analysis is the same as the one applied above to 
x(T) and £(T), for which the corrections to scaling are 
quite small. Namely, each critically-divergent observable 
F{j3) is represented by F opt (f3) of Eq. © over the whole 
range of (3 to a good approximation. The input con- 
sists of F*(/3), a confluent correction term and a very 
limited numbers of non-critical correction terms derived 
from HTSE. 

Together these results can be taken as validating the 
"extended scaling" approach. The approach could be 
systematically implemented in numerical work so as to 
improve yet further the accuracy of critical parameter 
sets derived for standard systems, possibly incorporating 
where necessary further higher order correction terms. 

Perhaps a more fruitful application would concern the 
analyses of numerical data in more complex systems, 
where the present accuracy of the critical parameter sets 
is much poorer. For instance, it has been pointed out that 
for the analysis of data on spin glasses with symmetric 
interaction distributions j3 should be replaced by (3 2 in 
all expressions^^ as all terms in the HTSE in these spin 
glasses are strictly even in fi. The extended scaling proto- 
col allowing for this and with appropriate 4>f(P) normal- 
ization factors has indeed been shown to significantly im- 
prove the consistency of critical exponent values derived 
from numerical simulations on Ising spin glasse a 6 ' 24 . 
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